Boundary lubrication with a glassy interface 
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Recently introduced constitutive equations for the rheology of dense, disordered materials are 
investigated in the context of stick-slip experiments in boundary lubrication. The model is based 
on a generalization of the shear transformation zone (STZ) theory, in which plastic deformation 
is represented by a population of mesoscopic regions which may undergo non affine deformations 
in response to stress. The generalization we study phenomenologically incorporates the effects of 
aging and glassy relaxation. Under experimental conditions associated with typical transitions from 
stick-slip to steady sliding and stop start tests, these effects can be dominant, although the full STZ 
description is necessary to account for more complex, chaotic transitions. 
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I. INTRODUCTION 

Advances in developing nanometer scale technologies 
and devices are intrinsically coupled to fundamental 
progress in scientific understanding of properties of ma- 
t^a^j^^nd^ ^ ^ P]fj ^ 

friction, fracture, and plastic deformation at microscopic 
scales influences the operation and performance of en- 
gineered systems. On the other hand, investigation of 
clean, well characterized nanoscale systems is leading to 
new insights into these underlying physical phenomena 
which govern systems under stress. 

Increases in scientific computing capacity along with 
the development of new experimental techniques, has re- 
cently created opportunities for progress in both theory 
and measurements. The surface force apparatus (SFA) 
was originally designed to study solvation forces induced 
by a liquid confined between parallel surfaces 0,0. Sub- 
sequently it was adapted to measure shear forces |j| • The 
SFA allows precise measurements on a microscopically 
thin layer of lubricant, separating atomically smooth 
(typically mica) surfaces. Friction and/or adhesion as- 
sociated with a single asperity contact can be precisely 
measured. We focus here on friction and the associated 
stick-slip instabilities. 

In many practical instances friction involves rough ma- 
terials. However, for rough surfaces it can be difficult to 
identify underlying mechanisms associated with complex 
phenomena (e.g. irregular dynamics, and bifurcations). 
For rough materials individual asperity dynamics are a 
priori combined with any collective phenomena which 
may be associated with the population of contacts, and 
the population of contacts itself is necessarily time de- 
pendent. Rough surface measurements are thus naturally 
complementary to investigations of plasticity and rheol- 
ogy of isolated, individual asperities. In addition, for sin- 
gle asperities, effects associated with interfacial materials 
(e.g. lubricants) are more easily isolated. Furthermore, 



there is growing evidence that the frictional properties of 
rough surfaces even at macroscopic scales is controlled by 
the plastic deformation of individual contacts. Indeed, 
experiments have recently been designed which isolate 
and measure the dynamics of individual asperities at a 
rough dry interface, subject to shear pjij] . 

In this paper we focus on dynamics of an individual 
lubricated asperity contact. We model the lubricant us- 
ing a set of constitutive equations which generalize the 
shear transformation zone (STZ) theory for amorphous, 
glassy materials. In addition to the STZ equations, the 
model incorporates the effects of glassy relaxation via 
the introduction of a state variable related to the inter- 
nal free-volume (the additional state variable may alter- 
natively be thought of as an out-of-equilibrium effective 
temperature) . The coupled STZ and free- volume dynam- 
ics were introduced previously, and shown to capture a 
range of experimental phenomena in glassy and granular 
materials |35l HH ■ Here the constitutive equations model 
internal states of the lubricant. We perform a series of 
analytical and numerical calculations which mimic typ- 
ical SFA experiments fl5t Eil l3fj| . and investigate the 
stationary states, the bifurcation diagram of the transi- 
tion between stick-slip and steady sliding, the nature of 
this transition (i.e. super- or sub-critical), the emergence 
of chaos, and aging of the yield stress in stop-start tests. 
At this stage, we primarily map out qualitative behaviors 
of the model. The compelling correspondence to exist- 
ing experiments sets the stage for future, more detailed, 
quantitative comparison with data. 

The remainder of this paper is organized as follows. 
Section II provides a brief overview of friction and bound- 
ary lubrication. Section III describes free-volume consti- 
tutive equations and summarizes elements of STZ theory. 
Section IV contains the numerical and analytical results 
in scenarios representative of SFA experimental studies. 
Finally, we conclude in Section V with a discussion of our 
results, comparisons with experiments, and directions for 
continuing research. 
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II. BACKGROUND 

In this section we provide a brief overview of boundary 
lubrication. Our emphasis is on novel material properties 
and the associated modeling challenges which arise for 
atomically thin, confined liquid films. Note that, even 
for interfaces with relatively simple features, we still lack 
precise quantitative, predictive models for friction. A 
more complete overview of recent results on friction and 
lubrication can be found in |12| and references therein. 

In boundary lubrication a molecularly thin film of ma- 
terial is confined between two parallel surfaces. The rela- 
tive motion of the surfaces is mediated by the plastic de- 
formation of the interfacial material lubricating the con- 
tact. Molecularly thin films display specific properties 
which differ from the viscous behavior of bulk materials. 
Among the most noteworthy of these is the development 
of a yield stress for thin films at sufficiently low temper- 
atures. 

The changes of material properties under confine- 
ment can largely be attributed to a liquid-solid transi- 
tion [5|,|8(. In some cases, numerical and experimental 
evidence s uggest s lay ering occurs in the interfacial mate- 
rial jlllllSlilEllIllil and the situation is similar to 
a liquid-crystal transition. That is, the material orders, 
and its deformation is expected to result from propaga- 
tion of dislocations, or layer-over-layer sliding. In other 
cases, no evidence of ordering of the interfacial mate- 
rial is observed, and the liquid-solid transition induced 
by confinement enters the very large class of structural 
glass transitions ITU E2I l32j . In this case the material 
remains amorphous (liquid-like), but displays solids- like 
properties. The emergence of a yield stress is accompa- 
nied by power-law viscosities, as well as signs of glass y ag - 
ing [28ll3oL l3^| and anomalous response spectra |22ll24|. 
For some lubricants, either a glassy or layering transi- 
tion can be observed, depending on features such as tem- 
perature [22, holding time, surface roughness , and 
commensurability of the surface and the film 23] . Con- 
sequences of solid-like ordering within the film include 
development of a static yield stress, and stick-slip insta- 
bilities 0E3. 

Boundary lubrication and the SFA experiments pro- 
vide special opportunities for theory, because the inter- 
face is well characterized and precisely controlled, yet the 
system is large enough to display phenomena which also 
arise macroscopically. To date, models have primarily 
emphasized effects associated with ordering and interac- 
tions between the crystalline substrate and the lubricant. 
These aim to describe effects of layering and surface in- 
duced order and involve (i) simple, reduced models of 
non-interacting particles in an effective periodic potential 
induced by the surfaces ^2j^ , (ii) Ginzburg-Landau func- 
tionals which account for heterogeneous ordering or 
(iii) motion of adsorbate layers in the periodic potential 



associated with a regular surface 0, 0, EH] . 

Alternative approaches focus on internal properties of 
the lubricant, including contrasts between liquid and 
solid-like properties and glassy behavior. Phenomeno- 
logical rate and state friction laws have been intro- 
duced [2(1 H(|, in which the friction depends on the 
instantaneous slip rate and a state variable. The state 
variable models the collective dependence of friction on 
the internal degrees of freedom of the lubricant. This 
approach assumes the fluctuations are sufficiently self- 
averaging that microscopic degrees of freedom in the 
boundary layer can be ignored. This simplifying assump- 
tion was inspired by a large body of work in dr y fr iction 
where rate-and-state formulations EH EH DH E3 have 
shown to be useful to account for exp erimental data in- 
cluding stick-slip instabilities EH El EH EH E3 . State 
variables can be motivated by experimental observations 
or molecular dynamics. For dry friction the state vari- 
able is related to the average lifetime of individual con- 
tacts, whereas in boundary lubrication the state variable 
is loosely connected to the degree of internal order in 
the lubricant. However, the friction laws these underly- 
ing mechanisms inspire are based more on macroscopic, 
thermodynamic-like criteria, than on the underlying mi- 
croscopic physics at the interface. 

Our model aims to provide a microscopically motivated 
description of the macroscopic forces which arise when 
the amorphous, interfacial lubricant is subject to shear. 
It relies on the assumption that the deformation of amor- 
phous materials is controlled primarily by excluded vol- 
ume effects, which dominate over fine details of molecular 
interactions. This is an old idea, advocated by Struik in 
the 70's EH], and supported by striking similarities be- 
tween very different amorphous systems. As a result, a 
relatively simple account of visco-plasticity is expected 
to hold for wide classes of amorphous systems. In the 
conclusion, we discuss how the equations we study here 
may account for the behavior of not on ly lu bricants, but 
also sheared granular materials EH EH EH E3 ■ 

For the theorist, boundary lubrication has some spe- 
cial, simplifying features which in many respects make 
it an ideal template to study plasticity of amorphous 
materials. The interfacial lubricant layer is sufficiently 
thin that certain bulk phenomena, such as strain local- 
ization, appear to be avoided. As a consequence bound- 
ary lubrication is, in fact, one of the few experimental 
setups where homogeneous constitutive equations can be 
directly tested. In contrast, for glassy bulk materials de- 
formation organizes in shear bands, usually a few par- 
ticle diameters thick, and the strain-rate measured in 
an experiment averages over a non-uniform field. The 
origin of strain localization is poorly understood, but 
clearly a more complete understanding of simple homo- 
geneous flows must be established first. No shear band- 
ing is expected to occur in nano-confined films because 
the deformation is already confined at scales equal to or 
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smaller than that which would be expected for the shear 
band. The SFA experimental setup resembles closely the 
sheared strip used in recent numerical studies of relax- 
ation in glasses [H . 

There are, however, several persistent challenges as- 
sociated with the SFA. Firstly, an astonishingly wide 
range of phenomena have been observed in this system 
for different lubricants under different conditions. The 
range of behaviors remains a puzzle which is difficult to 
piece together in the absence of a systematic theory that 
clearly captures the most basic observations. The vari- 
ability of experimental observations might be attributed 
to inherent difficulties in the preparation of samples of 
any glassy material, due to the effects of aging. Cer- 
tainly capturing the full spectrum of properties is a long 
term target for theoretical models. Secondly, it is dif- 
ficult to increase the stiffness of the apparatus beyond 
~ 3500 N/m |5^]- A salient feature of boundary lubri- 
cation, which directly results from this finite stiffness, 
is the emergence of stick-slip instabilities at low veloc- 
ity H Q El El EE Hi In the most dra- 
matic cases, the transition to stick-slip is accompanied by 
chaotic behavior [3(| . The main consequence of the stick- 
slip instability (and of the finite scanning length of the 
experimental apparatus) , is that the SFA can provide sta- 
tionary data only for limited ranges of strain-rates. Thus 
for a wide range of parameters, information about the in- 
terfacial material can be gathered only through the obser- 
vation of instabilities and transient dynamics. Therefore, 
it is essential to treat this finite stiffness explicitly in any 
theoretical approach. Ultimately, rather than hampering 
our understanding of structural glasses, stick-slip insta- 
bilities and transient dynamics provide us much richer 
data than simple stationary states. These instabilities 
may help clarify important, general issues associated with 
the out-of-equilibrium properties of structural glasses. 

The model we study in this paper is an intermedi- 
ate statistical theory recently proposed jH HI EEJ for 
sheared structural glasses. Here "intermediate" refers 
to the fact our model falls between the thermodynamic- 
like, phcnomenological rate and state descriptions which 
have been proposed for friction at dry and lubricated in- 
terfaces, and an atomistic statistical mechanical descrip- 
tion of the lubricant, which takes into account detailed 
microscopic, molecular interactions. The model begins 
with shear transformation zone (STZ) theory |56t l57| . 
which was inspired by molecular dynamics simulations 
of fracture in amorphous materials, and which pro- 
vides a microstructural description of shear-induced re- 
arrangements [H STZ theory accounts for 
the emergence of a yield stress in amorphous materi- 
als through introduction of state variables characteriz- 
ing the anisotropy of structural arrangements. Drawing 
on earlie r appro aches to describing creep in metallic al- 
loys mUfliiilii, STZ theor y models local shear 

rearrangements as activated processes, controlled by lo- 



cal density fluctuations. This leads to the introduction 
of free-volume activation factors, in the spirit of early 
theories of the glass transition [6J, |63, |6y, |63, |68| . 

In [35| and subsequent works [H], it was noted that 
previous approaches treat free- volume as a fixed param- 
eter, although it clearly varies as the material dilates or 
contracts. In granular materials, for example, density re- 
laxations have been observed and precisely characterized 
experimentally | 6 3 | , inspiring several models of slow re- 
laxations |69l l70l l7ll| . Dilatancy in granular material is 
also involved in the definition of frictional properties and 
stick-slip transitions [j^. Slow relaxation of volume or 
enthalpy are ubiquitous features accompanying the glass 
transition ^1] ■ We expect more explicit modeling of the 
time-dependence of density fluctuations may capture the 
emergence of a wide range of glassy properties near a jam- 
ming transition. The constitutive equations introduced 
in [23 were thus developed to address the question: How 
much glassy phenomenology can be captured by the sim- 
plest account of free- volume dynamics, coupled to the 
dynamics of shear transformation zones? 

Previously it was shown that free-volume dynamics 
suffices to characterize aging and power law rheolo- 
gies, while the dynamics of shear transformation zones 
are required to account for the emergence of a yield 
stress 36]. In transient regimes, both processes may in- 
teract and contribute to rheological properties. Further- 
more, chaotic behavior has been observed in the SFA 
close to the stick-slip instability [13], and it is known in 
over-damped frictional equations, that two or more state 
variables are necessary to understand the occurrence of 
chaos. Here, we did not invoke an ad hoc theory with 
multiple state variables [j^, by were led to it by the 
underlying physical mechanisms already associated with 
other, relevant experimental observations. Moreover, as 
chaos is difficult to characterize experimentally, signifi- 
cant insight is provided when observations are supple- 
mented with theoretical models, to guide measurements 
and analysis. 

Finally, we note that it is not a priori necessary to 
identify the additional state variable as a free-volume. 
It was previously noted by Falk and Langer that free- 
volume was related to the notion of Edwards' tempera- 
ture jH El El El As noted later [H, the essential 
feature that free-volume dynamics captures is the exis- 
tence of an intensive quantity which measures internal 
disorder, and which evolves as the system orders or is 
driven away from equilibrium. Alternative approaches 
characterize the internal state in terms of an effective 
temperature |55|. The similarity between free- volume 
and a temperature results from the fact that, in hard- 
sphere materials at constant pressure, the dominant con- 
tribution to the energy of a subset of molecules is en- 
thalpy. In this case, local energy and density fluctuations 
are directly correlated. Note finally, that the concept of 
Edwards' temperature is likely to be related to effective 
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temperatures arising in weak versions of the fluctuation- 
dissipation theorem EH Ell- I n this paper, we 
use free- volume terminology, and refer the reader to [5{| 
where the relation between free- volume and effective tem- 
perature is discussed in more detail. 



variables and/or nonlinearities, and to illustrate phenom- 
ena which require a larger number of state variables to 
be resolved. 

(ii) Elements of STZ theory 



III. CONSTITUTIVE EQUATIONS 

Our presentation of the model is broken down into four 
steps: (i) preliminaries associated with the SFA, explic- 
itly accounting for the finite stiffness, (ii) STZ theory, 
(iii) free volume equations, and (iv) rescaling to obtain 
dimensionless equations. The equations in sections (ii- 
iv) have been presented elsewhere, but we include their 
derivation for the sake of completeness. We refer the 
reader to j5|| and j3|| |55| for more detailed discussions 
of the underlying assumptions. 

(i) Surface Force Apparatus 

A primitive model of the SFA consists of a single slider 
block, pulled along a surface by a spring of stiffness k. 
The opposite end of the spring advances at a prescribed 
velocity V. Letting x denote the displacement of the 
block relative to the stationary surface, the spring ex- 
erts a force F = k{Vt — x), where t is time measured 
from some initial time t = when x — . The block 
is subject to both the pulling force from the spring and 
frictional resistance at the surface. Modeling the SFA 
as a single slider assumes that the sliding surface is suf- 
ficiently small, rigid, and uniform and that friction at 
the contact is sufficiently self-averaging that slip occurs 
uniformly across the interface. 

Assuming the thickness h of the interface remains con- 
stant, the motion of the slider is related to the rate 
of shear deformation 7 of the interfacial material by, 
x = 2 h 7. If the area S of the contact is constant, the 
shear stress exerted by the slider on the interfacial ma- 
terial can be written, a = F/S. Furthermore, for the ex- 
periments we consider, the friction is sufficiently strong 
that the motion is overdamped, and inertial forces associ- 
ated with the nonzero mass of the slider can be neglected. 
This leads to an equation of motion for the stress: 

<r = /i(e-7) (1) 

with, n = hk/S, and e = V/(2h). The motion x(t) of 
the slider follows from the solution to (1), which depends 
on how the strain-rate 7 is related to shear stress a and 
the internal state variables. 

We next define the constitutive equations which couple 
free- volume and STZ dynamics. Later we restrict our dis- 
cussion to subsets and (linearized) simplifications of these 
equations to isolate effects associated with distinct state 



STZ theory is based on the idea that the macroscopic 
deformation of an amorphous material results from local- 
ized rearrangements involvin g coope rative molecular mo- 
tion at mesoscopic scales [8ll Hj, |83| . The loci of such re- 
arrangements are called shear transformation zones, and 
the internal state of the system is characterized by their 
number density. In its simplest form, STZ theory involves 
only two types of zones (labeled "+" and "-" , with num- 
ber density n + and rt_, respectively) oriented along the 
principal axes of the shear stress. Zones of each types 
transform into one another during an elementary shear. 
The average strain rate 7 is given by 

7 = M{R+n+- R-n-) . (2) 

Here 7 averages over the populations of zones, n±, that 
reorient with probabilities R±, respectively. 

The important insight Falk and Langer contributed 
to previous STZ theories was to treat the populations 
densities n± as state variables, and propose equations 
of motion for the populations. These take the following 
form [Hil: 

n± = — R T n T + R± n± + cry (A c — A a n±) . (3) 

The first two terms on the right hand side account for ex- 
changes between the populations of STZ's due to meso- 
scopic rearrangements, while the last term introduces a 
coupling with the mean flow, through creation of STZ's 
at rate A c , and annihilation at rate A a - The equations 
describe how shear deformations induce small displace- 
ments of the molecules, hence creating and destroying 
shear transformation zones. In this framework, the emer- 
gence of a yield stress in amorphous solids at low temper- 
ature is associated with the mobilization of zones when 
stress is applied. 

(iii) Free-volume activation 

In the original formulation of STZ theory, the rates 
were estimated to be nonlinear functions of stress. 
Derivation of the rates were based on free-volume ac- 
tivation, as developed by Cohen, Turnbull and cowork- 
ers to understand the phenomenology of the glass tran- 
sition E3> Ell ■ Specifically, the rates R± were esti- 
mated to be of the form exp(— vo(±a)/vf), where vq is 
a stress-dependent activation volume, and Vf is a ma- 
terial dependent constant. The detailed formulation of 
transformation rates R± is not essential to capture the 
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STZ mechanism for jamming. Instead, a first order ap- 
proximation for the stress dependence of the rates R± is 
sufficient [H, E3l 

However, for hard-sphere systems, the free- volume v/ 
is directly related to the density. There is no reason to 
believe that it should take a fixed value as a function of 
pressure and temperature. On the contrary, it is a dy- 
namical quantity which evolves as the material dilates 
or contracts. This observation naturally leads to dy- 
namical equations for the free- volume, written in analogy 
with equations of motion for the populations n± |35ll36|. 
While free-volume accounts for the existence of disorder 
in molecular packing, the difference in STZ densities de- 
scribe its anisotropy. 

Activation factors depends both on free-volume and 
stress fluctuations. Assuming these effects are uncorre- 
cted, we express the rates as: 



R±(c, v f) = Ro ex P 



vo 

v f 



exp 



±- 



(4) 



An elementary shear rearrangement takes a "+" oriented 
zone to a "-" oriented zone, or vice versa, and occurs if 
sufficient free- volume (large than vq) is available, and if 
the virtual work of shear forces promotes the transition 
in the ± — + =F direction. The variable p, is here a scale 
of forces (not to be confused with the elastic modulus fi) 
which may depend on temperature, and governs stress 
activation factors. Note that the introduction of back- 
ward and forward jumps is an old idea, already present 
in Eyring's theory of viscous liquids l84l or in Spaepen's 
approach to creep in metallic glasses |5j| ■ 
The equation of motion for Vf is given by 



Vf 



-Ri exp 



Vf 



+ A v a A f 



(5) 



Here free- volume dynamics involve two competing mech- 
anisms: (i) activated elementary compaction which in- 
crease the density, and (ii) the transfer of macroscopic 
work into enthalpy, which dilates the material. The pa- 
rameter A v specifies how efficiently the work of external 
forces is used in dilatancy. The parameter R\ is an up- 
date frequency, which should be of the same order as Ro; 
the activation volume v\ may differ from vq , because the 
two elementary rearrangements (shear and compaction) 
involve different relative motion of the molecules, hence 
different reactional pathways. The ratio k = vi/vq is 
an essential parameter of the theory, and is expected to 
depend on the shape of the molecules. 

The equations described here correspond to the low 
temperature limit of a more general set of equations for 
the dynamics of a disorder temperature |55| . We will 
restrict our current discussion to this free-volume for- 
mulation. As we will show, it captures a wide range of 
phenomena observed experimentally. 



(iv) Rescaling and change of variables 

For the constitutive equations defined above (equa- 
tions 01 0), it is convenient to introduce reduced 
variables |56ll57l| . 



A 



A 



and 



Vf 

x = — 

(6) 

along with the rescaled parameters, = 2A C /A a , £q = 
AoA c /A a , Aio = 1/(A A C ), E = 2e i? , and E x = 
Ri/vq. This change of variables leads to the following 
set of equations: 



7 
A 

A 

x 



Eq exp 



A sinh ( — 



a j 



1 fl-^A 



Mo eo 
—Ei exp 



Ik) 

(I -A) 



acr7 



A cosh ( -=)p) 
(8) 
(9) 
(10) 



Equations 01, 01, and 01 - free- volume x m {El be- 
ing held constant - are almost identical to the ori gina l 
formulation of STZ theory by Falk and Langer [561 l57j . 
The variable A accounts for the total density of STZ's. 
The steady state value A = 1 is a stable fixed point. In 
STZ theory, A differs from 1 only in the initial transient, 
and its initial value is expected to depend on the type of 
annealing performed during the creation of a glass from 
a high temperature liquid. Here we are not concerned 
with transient features associated with these annealing- 
dependent initial values of A, and will assume A = 1 
throughout. This eliminates equation 01. The variable 
A, is a normalized difference "+" and "— " STZ densities. 
It accounts for the anisotropy of the molecular structure. 
For a fixed applied stress a (and fixed free- volume, x)j 
equations and were shown previously to account 
for a transition between an elastic regime (jamming) and 
a visco plastic regime (flowing) . The transition occurs at 
a yield stress, a y , satisfying 



tanh f ^ ) = £E 



(11) 



The variable x is a normalized free- volume, which ac- 
counts for the existence of disorder in the molecular struc- 
ture. Its dynamics, determined by equation (|f 0(1 . couples 
with equations 01 and only as far as x appears in 
the prefactor of equation 0, setting the time scale of 
elementary shear events. 

At a fixed applied stress a, equation 01, 01 and (fTTTf) 
account for plastic deformation resulting from the cou- 
pled dynamics of A and x- When the system is driven at 
a constant shear rate as in the SEA, these equations must 
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be supplemented with equation QJ, which accounts for 
coupling of material deformation with a compliant driv- 
ing apparatus. 

Discussion 



the transitions to stick-slip at low velocities. In the 
last part of our work, we focus on the existence of (iv) 
chaotic regimes of stick-slip. In order to observe chaos, 
we need the full (three-dimensional) set of non-linear 
equations (|7I1UI) . 



Most of the results presented here are consequences of 
free- volume dynamics IjlOjl coupled to equation Q, the 
STZ variables A (and A) being held constant (A = and 
A = 1). In this case, the equations J7J and (|10fl reduce 
to: 



7 



X 





r 




<7 


Eq exp 




sinh 






x_ 




A 



-Ei exp 



+ a a 7 



(12) 
(13) 



These equations can also be linearized for small stresses, 
which enables us to eliminate the parameter p, (taken to 
unity) which enters the exponential activation factors to 
fix a scale of stresses. The linearized (in a) equations are 



7 = E exp 
X = -Ei exp 



+ a a 7 



(14) 
(15) 



We refer to this as the stress-linear approximation. It 
captures most of the phenomenology accompanying the 
stick-slip instability, and allows for interesting analytical 
calculations. 

Throughout this paper, we will use these three sets of 
equations- the full, coupled system of equations includ- 
ing STZ and free-volume effects, the nonlinear equations 
describing free volume dynamics only and the linearized 
version of the free-volume equations- in order to clarify 
the consequences of our assumptions and the role of the 
different state variables. The questions raised are: What 
behavior is already captured by free- volume dynamics, in 
the stress-linear version (equations (|1411 5|) ^ ? What is the 
importance of activation factors, and the non-linear de- 
pendence of strain-rate versus stress (equations H12I13I) )? 
What are the expected consequences of the interaction 
between several internal state variables, as modeled by 
equations Q, JSJ and, (fT0f> ? 

IV. COMPARISON WITH EXPERIMENTS 

In this section, we compare the qualitative behavior 
of our model with several important experimental ob- 
servations. Most of our investigation focuses on the 
equations governing free- volume dynamics, either in their 
non-linear form (|12I13|1 or in their linear form l|14llfljl . 
We study (i) stress versus strain rate relations in steady 
sliding, (ii) transient dynamics upon start-up, and (iii) 



(i) Steady sliding 

The first step in characterizing the behavior of sheared 
materials is, of course, the steady state relation between 
stress and strain rate. In boundary lubrication and single 
asperity friction, it has been observed that, at low veloci- 
ties, the friction force is weakly velocity weakening. This 
weak dependence is reminiscent of the original observa- 
tions by Dieterich 0, in the case of dry surfaces, 
where a logarithmic dependence can be assumed with a 
reasonable degree of confidence. In boundary lubrication 
and single asperity friction there is insufficient data to 
conclusively determine whether the steady state relation 
between force and velocity is logarithmic, power law, or 
a combination of the two [2^, Below we calculate 

the stress vs. strain rate relationship in order of increas- 
ing model complexity, beginning with the linearized free- 
volume equations, followed by the nonlinear free volume 
equations, and finally for the complete model, including 
STZ's. 



Stress- 



equations 



In the stress-linear version of our model (equations i|14|) 
and flU ) the ratio between stress and strain rate deter- 
mines a viscosity, 



77 = exp(l/x)/£o 



(16) 



which is a simple function of the free-volume x- 

Initially, we take the limit of infinite stiffness for the 
experimental apparatus, so that, from equation 7 = 
e. ^From equations 114J1 and l|15(l . the dynamics of x 
reduces to: 



X 



-Ei exp 



E 



exp 



(17) 



The quantity x admits a steady state value only if the 
shear rate is not too large. For high shear rates, 



Eq Ei 



(18) 



this equation becomes unstable, x is positive at all times, 
and x diverges. This instability is not related to stick- 
slip, since it occurs for any value of the stiffness. It indi- 
cates that at large shear rates, the material cannot dis- 
sipate the work of external forces, and is driven towards 
a highly disordered state. Situations when e > e* are by 



7 



definition transient: the material does not reach a steady 
state. Instead x diverges. 

Interpreting \ as a free-volume, this divergence may 
correspond to an opening of the interface, and possible 
loss of contact between the surfaces. In typical experi- 
ments, if such high shear rates are applied, the divergence 
of x would be a priori limited by the scanning length of 
the experimental device. It is intriguing, however, that 
opening modes at large shear rates are observed experi- 
mentally |53| . They have not been characterized in the 
literature because they are usually viewed as experimen- 
tal technicalities which limit the range of accessible driv- 
ing velocities. Instead, experiments have been restricted 
to situations in which such effects do not occur, specifi- 
cally low driving shear rates, e < e*, or high shear rates 
on limited scanning lengths. In the latter case, the appar- 
ent viscosity, rj = exp(l/x)/£'o ~ V-Eo is approximately 
independent of x, and it may appear, from a measure- 
ment of stress alone, that the system is stationary al- 
though the internal dynamics may not have reached a 
steady state. 

The existence of this divergence at a critical driving 
shear rate e* arises physically from the fact that the tran- 
sition probabilities are bounded. The factor exp(l/x) 
approaches unity at high %, and becomes decreasingly 
sensitive to changes in the intensive variable x when x is 
much larger than 1. At this point, free- volume dynam- 
ics is only weakly coupled to the other equations. Of 
course, if a high shear rate e > e* is applied steadily, this 
divergence ultimately leads (in the t — > 00 limit) the sys- 
tem towards a highly disordered, "fluidized" , state. Our 
equations are not designed to describe this limit, and in- 
stead should start to break down when x reaches some 
large values, say Xf- We expect Xf >> 1 (here the value 
1 corresponds to the activation free-volume for an ele- 
mentary shear transformation, which is much less than 
one molecular volume per molecule, whereas Xf is of the 
order of one additional molecular volume of free space 
per molecule). We do not attempt to account for the 
late stages of this divergence, and only changes of mate- 
rial behavior around x ~ 1 are of interest to us. For this 
reason, and to simplify the discussion, we can safely take 
Xf to infinity. 

We could avoid explicit reference to such a divergence 
as several mechanisms could be invoked to explain the 
saturation of free- volume at high velocity resulting in the 
emergence of another fixed point. Possible mechanisms 
include: weak density dependencies of various parame- 
ters, like Eq, Ei, or a, or higher order terms in the free- 
volume dynamics (equation (fTUl) . l(T5|) . or (|T3|l ). Here 
we study free-volume dynamics in the simplest form for 
several reasons. Firstly, several mechanisms could be in- 
voked to account for the saturation of free- volume at high 
velocity, and it is unclear at this stage which would be 
dominant. Thus we prefer not to differentiate between 
them. Secondly, regardless of the mechanism, the ap- 



proach to any high free-volume fixed point would, in its 
early stages, closely resemble the divergence we discuss 
here. Given the limited span of experiments, we expect 
our model is relevant to observations which are currently 
available. 

For shear rates smaller than e*, free- volume admits a 
steady state value, and this leads to a relation between 
stress and strain rate, of the form: 



e 



EqEi 



(19) 



The constitutive equations thus account for power-law 
viscosity, with an exponent 



(20) 



1 



which is directly related to the ratio k = vi/vq of ac- 
tivation volumes. For k > 1 the stress is an increasing 
function of the shear rate. The material is shear thicken- 
ing. For k < 1, the stress is a decreasing function of the 
shear rate, the material is shear thinning. We will see 
that, in the later case, the system exhibits a transition 
to stick-slip at low velocities. 



Non-linear free-volume equations 

Next we calculate the stress vs. strain rate for the non- 
linear free volume equations Ijl2(l and (|13J) . In this case, 
we obtain a generalization of Eq. (|19fl which is valid be- 
yond the range of small stresses, where the linear stress 
approximation is expected to hold. The resulting rela- 
tionship is given by: 



/ olEq 



V E 1 



a ( sinh [ — 



(21) 



For small a linearization of the hyperbolic sine, leads to 
the power law rheology described above (Eqs. I|19H20|) '). 
For large stresses, the right hand side is dominated by 
the exponential growth in the hyperbolic sine, which re- 
sults in a logarithmic dependence of the stress a on strain 
rate e. This expression is compatible with experimental 
observations, although as already mentioned, the range 
of available observations is sufficiently limited to render 
detailed functional fits inconclusive [3, 153, • 

As in the linearized case, Eq. I|21(l is accompanied by 
a condition imposed on the strain rate e which must be 
satisfied for there to be a self-consistent, steady state 
value of x- Violation of the condition (e > e*) is associ- 
ated with diverging free volume and fluidization. For the 
nonlinear free volume equations the constraint is more 
complicated than in the linearized case, and is described 
by e/(Eo sinh(cr//2)) = exp(— l/x) < 1- A stationary 
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value of x exists if and only if the inequality is satisfied. 
Equivalently, a solution exists iff: 



( aE . I fj 
■ a sinn 



V E 1 



< 1 



(22) 



The case of equality in Eq. i12_?ll can be used to define a 
critical value a* : 



aE .-,(<?* 
a sinn — 



I 



(23) 



Then, validity of Eq. (|22|l depends on the value of a rel- 
ative to er*. If k < 1, er is a decreasing function of e, 
and the condition 11' I'D is met for er > er*, or equivalently, 
e < e* = e(c*). If k > 1, er is an increasing function of 
e, and the condition l|22|) is met for er < er*, or e < e*. 
In both cases, as with the linear equations, there is a 
limiting driving strain rate, e* above which the material 
fluidizes, and (in our model) the high- velocity fixed point 
is sent to infinity. For e > e*, as the quantity x diverges, 
the relation between stress and strain rate reduces to 



(24) 



e = Eq sinh — 



Stress appears to be increasing for large values of the 
strain rate. 



(ii) Transient Dynamics Characterizing the 
Approach to Steady Sliding 

Next we describe transient effects associated with dis- 
crete jumps in the applied strain rate. Both the initial- 
ization of the system, as well as several common exper- 
imental procedures probing transients under controlled 
conditions, can be described in this general framework. 
We first consider the initial, waiting time dependent tran- 
sients associated with starting the system from rest. This 
is followed by an examination of the experimental proto- 
col, referred to as stop-start, or slide-hold-slide tests, in 
which the system is prepared in the steady state before 
the waiting time begins. For this analysis, we will focus 
on the linear and nonlinear free volume equations. Inclu- 
sion of the STZ terms does not quantitatively alter the 
results. 



Transient dynamics upon start-up 

In the absence of forcing, a = 0, the linear and non- 
linear free volume equations are equivalent. /From equa- 
tion (|I0|l (or (|I3[I . or H15|) ~). the free volume \ relaxes to 
0. The late stages of this relaxation (as the waiting time, 
t w , goes to oo) are logarithmic in time: 



Full Free Volume and STZ equations 

Finally, we explore the steady state stress vs. strain 
rate relationship using the complete set of equa- 
tions 0, (JBJ, ©, and {TO}. In this case, the steady 
state relation between stress and strain rate reads: 



Eq 



( aEo <J 



( sinh ( — ) cosh ( — 

V W a VM 



(25) 

The complete equations capture the physical phenomena 
of jamming (i.e. the absence of flow), which occurs for 
stresses smaller than a yield stress, <r y , which satisfies 
equation Ijlll) . Otherwise, the behavior is very similar to 
our previous results for the nonlinear free volume equa- 
tions. For non-vanishing shear rates, we can again define 
er* > a y as the solution of, 



aE 
Ex 



sinh ( 



fio cosh — 



= 1 (26) 



and e* = e(c*). Again, for n < 1, the stress vs. strain 
rate relation is decreasing, and for k > 1, it is increasing. 
Again, the range of strain rates where steady sliding is 
reached is bounded by e < e* . We see from equation l(2*B|) 
that the STZ dynamics do not significantly change the 
sliding properties as soon as the relevant values of the 
stress are large compared to a y . This situation is partic- 
ularly relevant for stick-slip instabilities, which occur for 
k < 1, in the low velocity, large stress regime. 



log^iU) 



(27) 



Simultaneously, as the system becomes increasingly com- 
pact, the effective "viscosity" (Eq. increases: 



V(tw) = ~et exp 

Eq 



"^(£lU 1/K 
E 



(28) 



The immediate consequence of this time dependent re- 
laxation is that the material displays an age dependent 
initial value of the stress when shear is subsequently ap- 
plied. 

This is illustrated in Figure ^ where we plot \ an d 
er as a function of time, obtained from the integration 
of equation coupled to HI2|) and (|f 3fl . Distinct peaks 
correspond to samples of different ages at the onset of the 
applied strain rate. In every case, the system is initialized 
at t = with a large value (x(t = 0) = 10) of free- volume. 
^From time t = to t = t w , free- volume relaxes in the 
absence of applied shear. At time t = t w , a strain rate 
e is suddenly applied, and held constant from that time 
on. For each curve, the value of t w corresponds to the 
time when stress begins to ramp up. 

The initial increase of the stress is due to the fact that 
X has relaxed to a relatively small value during the wait- 
ing time, resulting in a high effecting viscosity, which re- 
sists rapid deformation of the material. When the stress 
gets large, free- volume increases. The system dilates, and 
the shear deformation rate increases suddenly, resulting 
in the observed stress drop. 
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FIG. 1: Transient Dynamics upon start up from the in- 
tegration of equation coupled to fl I Lift and with 
£7o = -Bi = a = 1. From equation 118H . e* = 1, and the 
applied strain rate is e = 0.9. Results for k = 0.8 (solid lines) 
and k = 1.2 (dashed lines) are shown. No stress is applied 
during the initial density relaxation of the material from t = 
to t = t w , marked by the increase of stress. At t w the strain 
rate is suddenly applied and the ensuing dynamics of 1/x 
(top) and a (bottom) are displayed. Here smaller values of k 
(i.e. k — 0.8) leads to a larger value of the dynamical yield 
strength, which results from the fact that a smaller value of \ 
(more compact state) is reached is reached during the waiting 
period. 



Stiff apparatus 

To provide analytical estimates of the age dependence 
of the peak stress, we take the limit of an infinitely stiff 
apparatus, \i — > oo. When a fixed shear rate is imposed 
after some waiting time t w , the interfacial material resists 
with a force which is an increasing function of age t w . 
For the linearized equation i[14jl. the stress upon start-up 
reads, 



hj 



1/k 



(29) 



The response is viscous, with an instantaneous initial vis- 
cosity which increases as a power of waiting time. 

For the nonlinear free- volume equations i|12|) the power 
law dependence is replaced by a logarithm. For large t w , 



the force upon start-up is: 



OStart (tw) — — log (Ext w ) + \X log ( -J- 



(30) 



This logarithmic dependence on waiting time and shear 
rate at long times is consistent with experimental ob- 
servations, although it must be noted that experiments 
performed with the SEA [2^ . Isij l are rather inconclusive 
and that most significant evidence comes from the sin- 
gle asperity experiments by Bureau etal j34|- Here, the 
logarithmic dependence is a direct consequence of the 
assumption that transformation rates R± depend expo- 
nentially on stress (i.e. proportional to exp(±c//2)). 



Stop-start tests 

The characteristics of glassy materials depends sen- 
sitively on sample age and preparation methods. It is 
therefore important to focus on experimental protocols 
in which the initial state can be relatively well-defined. 
One common experimental convention which aims to con- 
trol the initial state involves starting from steady sliding 
motion. The initial state of the material is determined 
by the driving velocity. In velocity step experiments the 
drive velocity undergoes discrete changes from one value 
to another, and the transient response is monitored. This 
was one of the original protocols used to investigate the 
correspondence between rate and state laws and dry fric- 
tion experiments [39t Ecjl. a nd was recently investigated 
for lubricated contacts 28] . Another protocol, referred 
to as stop-start or slide-hold-slide experiments, involves 
preparing the system in a constant velocity steady state, 
then suddenly stopping the drive and letting the system 
relax for a time t w , and finally restarting at the initial ve- 
locity. This latter protocol has been studied extensively 
in boundary lubrication 0,0,113' an d f° r the plastic re- 
sponse of single asperities at the contact between rough 
surfaces [j^. 

This stop-start protocol grants direct access to the ag- 
ing process from a well-controlled initial state. Stop-start 
tests are depicted Figure based on numerical integra- 
tion of equations (|12(1 and l|13|l in the limit of an infinitely 
stiff apparatus. Here x an d cr are plotted as functions of 
time. The initial value of x corresponds to the steady 
state at a given shear rate e = 0.5. When the exter- 
nal drive haults (the "stop" phase) x relaxes to a smaller 
value, as the lubricant becomes more compact. The same 
shear rate in then suddenly reapplied (the "start" phase) 
after different waiting times. In this /1 — > 00 limit, 7 = e, 
and the stress immediately takes its peak value, where 
the position of the peak marks the waiting time for each 
curve. After start-up, free-volume increases, due to the 
transfer of the work of shear forces into enthalpy. This 
results in a decrease of the viscosity 77, accompanied by 
a decay of the stress with time. 
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FIG. 2: Numerical integration of equations 1121 and I13H 
during stop-start tests for shear rate 7 = e = 0.5. Parameters 
are Eq — Ej = a = fi = 1, and « = 0.8. For each curve, 
the applied shear haults for a time t w , after which it is sud- 
denly reapplied. The dynamics of \ (top) and a (bottom) are 
displayed. The solid lines indicate the relaxation of \ in the 
absence of shear, and the value of a upon start-up. Different 
line styles are used for different t m . 



FIG. 3: Peak stress, Act as a function of waiting time t w , 
from integration of equations 112^ and 1)13(1 . Parameters are 
Eq = Ei = a = n = 1, and k = 0.8, whence e* = 1. Dif- 
ferent shear rates have been used: 7 = 6 = 0.3 (solid line), 
0.9 (dotted line) and 0.99 (dashed line). Approaching e* = 1, 
a latency time interval becomes increasingly well defined, in- 
dicating the increasing importance of free-volume dynamics 
beyond the small-x domain, where log-time relaxation of free- 
volume holds. 



The solid line in Figure marks the envelope of all re- 
sponse peaks, and thus defines the peak stress as a func- 
tion of the resting time t w . Experimentally, for short 
waiting times, Yoshizawa et al. 15]. observed that for 
certain lubricants, the corresponding curve exhibits a 
well defined latency time. That is, there is a thresh- 
old in the waiting time, below which no increase in the 
stress is observed. For longer waiting times, Drummond 
et al. and Gourdon et al. j2j|,|3i| found that the difference 
Act between the peak value of the stress, and the steady 
state value (the so-called "stiction spike") increases as 
Act ~ log(^) for large t w . Both the short and long 
time behavior is reproduced by the nonlinear free vol- 
ume equations, as shown in Figure |21 for different values 
of the driving velocity. Note the latency time (roughly as- 
sociated with the rapid rise of Ac) becomes increasingly 
sharply defined at high drive velocities. The follows from 
the fact that at high velocities, approaching the limit- 
ing shear rate e*, the stationary value of free- volume is 
large, hence free-volume activation factors (which scale 
as exp(l/x)) are essentially constant. Therefore, the free 
volumthc free volume x nearly decouples from the rela- 



tion between stress and strain rate. If the external drive 
is stopped during a short time interval, \ relaxes, but the 
changes in x have little effect on the dynamic viscosity 77. 
Therefore, when the shear is applied again, the relation 
between stress and strain rate is still very close to steady 
state, and only a very faint peak is observed. 

Note finally that, due to the small spatial extent over 
which shear can take place in the SFA, it is likely that in 
some experiments the internal state of the material never 
reaches the steady state value. It might appear that a 
steady state has been reached when x is large enough to 
almost decouple from the stress vs. strain-rate relation. 
In this case, the latency time depends on the time it takes 
for x to relax to values of order unity, which will depend 
on the details of sample preparation. This subtle history 
dependence may explain larg e fluctuations observed in 
some experiments l30| 
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(iii) Instability and transition to stick-slip 



Here, we show that our constitutive equations not only 
account for the existence of stick-slip behavior at low 
drive velocity, but that they reproduce the shape of the 
stick-slip cycle with remarkable accuracy. They also ac- 
count for the existence of continuous and discontinuous 
transitions to stick-slip depending on the stiffness of the 
apparatus, and in some cases, chaos in the neighborhood 
of the transition. For this analysis, the nonlinear free vol- 
ume equations will be sufficient to characterize the tran- 
sition from stick-slip to steady sliding, and the nature 
of the bifurcation. The full set of coupled constitutive 
equations, including STZ, preserve the phase boundaries 
described by the free volume equations. However, STZ 
effects (particularly, which introduce another dimension 
to the dynamical system) are required to capture chaotic 
phenomena. 



Stick-slip motion: shape of the pulse 

At low drive velocities, and for a sufficiently compli- 
ant apparatus, rather than sliding at constant velocity, 
the interface exhibits stick-slip motion. We begin our 
analysis of this motion by illustrating a typical stick- 
slip cycle obtained from numerical integration of equa- 
tion JJJ coupled with the constitutive equations H12|l and 
1)13)1 . In Figure ^ we plot the time series of x an d < J 
during stick-slip. As in the experiments, stick-slip cy- 
cles appear qualitatively similar to stiction peaks. How- 
ever, stick-slip arises when the driving motion is con- 
stant, unlike stiction peaks in stop-start test, which rep- 
resent transient responses to time varying slide-hold-slidc 
drivers and are typically monitored in regime where the 
steady state motion corresponds to constant velocity slid- 
ing, Stick-slip arises due to an internal instability of the 
material at a given shear rate. Rather than sliding at 
that steady rate, the material alternates between "stick- 
ing" and "slipping." During the stick phase free- volume 
decreases to a value which is small enough to hinder the 
relative displacement of lubricated surfaces, during which 
time the material creeps at a rate which is too slow to 
keep up with the external drive. Consequently, the stress 
builds up. When it becomes large enough to trigger di- 
latancy, the slip phase begins, and x increases suddenly, 
and the stress is released during rapid sliding motion. 

In Figure [S] the same data is represented in a plot of 
a versus 7, which compares quite favorably with typical 
cycles observed by Drummond and Israelachvili |28j- 



Characterization of the Hopf Bifurcation 

With decreasing velocity and for a compliant appara- 
tus, our constitutive equations exhibit a Hopf bifurcation 
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FIG. 4: Stick-slip dynamics obtained from numerical inte- 
gration of equations 1)12(1 and for a fixed strain rate 
e — 0.2, and stiffness fj, = 0.2. Parameters are Eo — E\ — fl — 
a = 1, and k — 0.8. The initial value of the free-volume is 
X = 1. The regime of steady plastic deformation is unstable 
and leads to stick-slip motion. Fast relaxations of the stress 
result from sudden dilatancy of the material. 



separating steady sliding from stick-slip dynamics. The 
locus of bifurcation points defines a phase boundary in 
the \i (stiffness) vs. e (strain rate) plane. A systematic 
analysis of the emergence of stick-slip motion is easily 
performed for the stress-linear version of the constitu- 
tive equations. The details are given in Appendix I. The 
critical stiffness defining the Hopf bifurcation point for 
Eqs. (JTJ, and {T!| is given by 



A^hopf - -rr 



E x 1 



E (« + l) 2 \E E 1 



In 



Eo Ei 



(31) 

Solving for /1 as a function of e with all other parame- 
ters held fixed, defines the phase boundary in the (fl, e) 
plane, below which steady sliding becomes unstable (see 
Figure ©). Our curve is exhibits the same convexity as 
available experimental data pjol Eif . Additional experi- 
ments are needed to provide a quantitative experimental 
characterization of the curve. 

Experimental data also suggest that different types of 
transitions between stick-slip and steady sliding are pos- 
sible. The transition may be continuous (approaching 
the phase boundary the amplitude of stick slip motion de- 
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FIG. 5: Stick-slip cycle from the data of Figure 2] Arrows 
indicate the direction in which the stick-slip cycle is followed. 
The ramps on Figure ^] correspond here to the increase of a 
at vanishing j. The plastic strain rate suddenly increases at 
almost constant a before the friction force starts to decrease. 



creases continuously to zero as drive velocity is increased) 
or discontinuous (there is an abrupt change from finite 
amplitude stick-slip spikes to steady sliding as drive ve- 
locity is increased). Continuous transitions correspond to 
supercritical Hopf bifurcations, and discontinuous tran- 
sitions correspond to subcritical bifurcations. The later 
case is typically accompanied by hysteresis (i.e. coexis- 
tence of stick-slip and steady sliding, such that a steady 
decrease in the drive velocity results in a transition to 
stick slip at a lower strain rate than that associated with 
the transition to steady sliding associated with increas- 
ing velocity from the stick-slip ph ase) . This feature is 
observed in both lubricated friction |3fJ as well as dr y in - 
terfaces. in the latter case it is better characterized |45j . 
though the physical origin of the similarity (if indeed it 
persists upon more detailed experimental investigations) 
remains unclear. For both cases, to date analytical mod- 
els have failed to capture the existence of both super- and 
sub-critical Hopf bifurcations. 

The transition point separating super- and subcritical 
Hopf bifurcations on the phase boundary in the (//, e) 
plane depends on third order terms in a normal expan- 
sion of the dynamical system around the steady state. 
Therefore, this property offers a particularly stringent 



FIG. 6: Phase diagram for equations Q, 1141 and I15H . and 
ft < 1, in loge-log/^ plane. The dotted line indicates the limit 
velocity e* beyond which x diverges. The solid and dashed 
lines denote the curve /ihopf below which steady sliding motion 
is unstable; /^hopf vanishes at e = e*. The solid part of this 
line corresponds to points where the transition to stick-slip is 
continuous (super-critical) while the dashed line corresponds 
to points where the transition is discontinuous (sub-critical). 
When the transition is discontinuous, there is a zone of bista- 
bility, which lies above the dashed line. 



test of constitutive equations. We have performed the 
analysis of the system (JJJ, H12[l and 1)13(1 . Our analysis 
defines a point: 



Merit 



_g e -K+v'2-2K+K 2 
£b 1 - K 



(32) 



on the Hopf line, which is drawn on Figure |H1 Above, 
and on the left of this point, the Hopf bifurcation is su- 
percritical (continuous); below, and on the right of it, the 
bifurcation is sub-critical (discontinuous). 

The shape of the phase boundary as well as the types 
of transitions, and even their relative placement in the 
phase diagram are consistent with current experimen- 
tal observations. However, a complete and quantitative 
characterization of the phase diagram for boundary lubri- 
cation remains an open challenge experimentally. Over- 
coming obstacles associated with the limited range of 
stiffness which can be probed with existing techniques 
will enable key observations which can be compared with 
both the contrasting case of dry friction experiments and 
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the theoretical results presented here. 

Chaotic stick-slip 

Drummond and Israelachvili [2^] observed that in some 
cases stick-slip motion became erratic in the neighbor- 
hood of the transition. These authors identified the er- 
ratic motion as chaotic, based on the analysis of experi- 
mental time series. Their results were indeed suggestive 
of positive Lyapunov exponents. However, methods to 
identify chaos based on analysis of an individual exper- 
imental time series are necessarily approximate and in- 
conclusive, and alone are insufficient to guarantee that 
the erratic motion is chaos as opposed to some form of 
noise amplification. 

Given these uncertainties, models can play a special 
role in helping to differentiate between mechanisms which 
may lead to chaos vs. other modes of irregular motion. 
Here we study the full set of nonlinear constitutive equa- 
tions Q, H7I10|) . which include the STZ state variables. 
The dimension of this dynamical system is three, which 
is the minimum requirement for chaos. As shown below, 
over a restricted range of parameters, the model admits 
chaotic solutions, characterized by irregular stick-slip and 
positive Lyapunov exponents. Because these equations 
are deterministic, noise amplification is ruled out as a 
source of irregular motion in the model. 

In order to measure Lyapunov exponents we directly 
integrate the equations of motion of a Lyapunov vector 
u(t) in the tangent space (see e.g. [8q'l. The largest Lya- 
punov exponent A is estimated from the long-time be- 
havior of the norm of this vector: 

A = lim - log||u(i)|| . (33) 

t — >oo t 

For this numerical study, we used a fourth order Runge- 
Kutta algorithm with fixed and variable time steps. We 
checked that the results do not change qualitatively as 
a function of the precision of our numerical method, as 
long as the precision is sufficiently high. 

Figure shows typical traces in the chaotic regime. 
The three variables a, A, x, as we U as the value of 

\{t) = - t \og\\u{t)\\ , (34) 

are shown. This data illustrates that chaotic stick-slip 
can arise in this system. Chaos results from the interplay 
between the dynamics of the variable A, which charac- 
terizes the anisotropy of the molecular packing, and the 
variable x which characterizes the dilatancy of this pack- 
ing. 

Further investigation reveals that the range of param- 
eters where chaos (i.e. at least one positive Lyapunov 
exponent) is observed is restricted to a relatively com- 
pact area of the \i vs. e phase diagram, which is close 



X 




2500 3000 3500 4000 





i 




i 1 



2500 3000 3500 4000 



a 




2500 3000 3500 4000 




t 

FIG. 7: Chaotic motion: A, x i a alm A are plotted as a 
function of time t for equations Q, Q, JHJ, I1UI . for param- 
eters, Eq = 0.5, E\ = 1, eo = 1, oe = 2, fio = 0.5, ft = 0.5, fi = 
0. 14, k = 0.8 and e — 0.14. The asymptotic limit reached 
by X(t) is the largest Lyapunov exponent, which is clearly 
positive, indicating chaos. 



to, but not overlapping the transition from stick-slip to 
steady sliding, in a region near the critical point sep- 
arating sub- and super-critical Hopf bifurcations. This 
relatively restricted range is not surprising. Since chaos 
requires the dynamical system to be three-dimensional, 
it may disappear as soon as one variable is enslaved to 
another. In our case, this decoupling may occur either 
when the time scale of the dynamics of A and x are well 
separated from each other or when the absolute value of 
A becomes negligible in equation Q. 

To illustrate the chaotic domain, in Figure|H]we present 
a two-dimensional image illustrating the value of the 
maximum Lyapunov exponent as a function of the pa- 
rameters \x and e. In the color map, negative values of 
the Lyapunov exponent appear in blue, and describe the 
case of a stable fixed point in the dynamics. This corre- 
sponds to steady sliding motion at constant x- This is 
the solution between the Hopf bifurcation and the e = e* 
line (recall for e > e* , x diverges, and our theory breaks 
down). In the color map, green corresponds to a zero 
value of the Lyapunov exponent. The value A = char- 
acterizes a time translation invariant manifold. Green 
thus describes stick-slip regimes, which exhibit a peri- 
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V. CONCLUSION 
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FIG. 8: Two-dimensional plot of the asymptotic Lyapunov 
exponent in the (log(e), log(/i))-plane obtained from numer- 
ical integration of equations (1, 7-10). The color scale on 
the right gives the maximum value of the Lyapunov exponent 
A, as defined in equation (40). Negative values of the Lya- 
punov exponent appear in blue: they correspond to steady 
sliding. The vertical asymptote of the phase boundary as 
log(/i) — > — oo (i.e. roughly the right edge of the figure) cor- 
responds to e*. Periodic stick slip corresponds to green, cor- 
responding to vanishing of the Lyapunov exponent. Chaotic 
behavior is associated with positive values of the Lyapunov 
exponent, which appear here in red and yellow. Chaotic re- 
gions of stick-slip are separated by windows of periodic mo- 
tion. 



odic limit cycle. Green also appears as the maximum 
Lyapunov exponent for large values of the driving shear 
rate (e > e*). In this case, the asymptotic state of the 
system is not strictly speaking a fixed point, because x 
is evolving at all times. 

Chaotic motion is associated with positive Lyapunov 
exponents which appear as red or yellow in Figure [5] 
They occur within the stick-slip portion of the phase dia- 
gram, near the transition between super- and sub-critical 
bifurcations. As shown in the figure, in the model this 
behavior is clearly separated from the Hopf boundary by 
a narrow range of regular stick-slip. The close proxim- 
ity of the chaotic zone to the Hopf bifurcation may ex- 
plain the fact that experimentally it appears that when 
chaotic stick-slip is observed, it merges continuously into 
the Hopf transition. While our equations predict a reg- 
ular Hopf bifurcation to periodic stick-slip, with chaotic 
stick-slip resulting from secondary instabilities, the inter- 
mediate periodic regime occurs over a sufficiently narrow 
range that it may be difficult to differentiate experimen- 
tally between this case and a case in which the Hopf 
transition takes place directly via chaotic stick-slip. 



We have shown that a limited set of constitutive equa- 
tions capture at least qualitatively a wide variety of 
phenomena observed experimentally in single asperity, 
boundary lubricated friction. Most of the phenomena 
are captured by the dynamics of a single state variable 
identified here as free- volume. In steady sliding, the free- 
volume dynamics account for the existence of velocity- 
weakening or velocity-strenghtcning friction laws, involv- 
ing mixed power-laws and logarithmic dependences. In 
stop-start experiments, logarithmic increases in the peak 
stress with increasing hold time results from the relax- 
ation of free-volume, whereas for short hold times a la- 
tency time emerges from the non-linearity of the transi- 
tion rates. The dynamics of free- volume also accounts for 
the existence of a transition from steady sliding to stick- 
slip as drive velocity is decreased, and for the presence of 
both continuous and discontinuous transitions, depend- 
ing on the stiffness of the driving apparatus. Finally, 
the complete set of equations where free- volume couples 
to the dynamics of STZs, accounts for the emergence of 
chaotic behavior close to the stick-slip transition. 

Our present work is in sufficiently close qualitative 
agreement with experiments to warrant a further round 
of more quantitative experimental vs. theoretical com- 
parisons. Our work leads to clear predictions for the 
friction in steady sliding regimes, and the transition to 
stick-slip as a function of the compliance. Available ex- 
perimental data provides steady sliding friction over a 
limited range of velocities, and a collection of stick-slip 
cycles at low velocities for a given interfacial material, 
most often for a unique value of the compliance. Avail- 
able measurements of power-law, or logarithmic relation- 
ships between force and driving velocity are not sufficient 
to validate or invalidate a set of constitutive equations. In 
principle, we could attempt to directly fit cyclic stick-slip 
data, but in practice, this is extremely difficult, requir- 
ing additional assumptions about initial values of one or 
more internal variables. What information can we thus 
use to compare theory and experiments? Compliance is 
included in the model at a low cost in terms of the as- 
sumptions required, yet provides a large panel of predic- 
tions which magnify the sensitivities of the constitutive 
equations and the microscopic assumptions of the theory. 
The location of the transition to stick-slip, its dependence 
on compliance, and the type of transitions to stick-slip 
are important predictions of our model that present op- 
portunities for more detailed comparisons with experi- 
mental data. Drummond and Israelachvili [2^] have been 
able to study transition to stick-slip for a few values of 
the compliance over a limited range. We hope this initial, 
exploratory study will be followed by a more complete set 
measurements to compare with our predictions. 

Recent experimental results by Gourdon and Is- 
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raelachvili 32] have focused on the temperature and pres- 
sure dependence of the transition to stick-slip. While 
these variables are more accessible experimentally than 
compliance, they present greater challenges for theory, 
and cannot be addressed directly here because temper- 
ature and pressure do not appear explicitly in the con- 
stitutive laws we consider. However, a recent derivation 
of similar equations provides new theoretical insights 
for temperature and pressure dependence, suggesting fu- 
ture opportunities to extend our current analysis in this 
direction. 

A quantity which arose in several contexts through- 
out this paper was the critical strain rate e*, beyond 
which the value of x cannot reach steady state. In our 
equations, above e*, x diverges with time, because the 
internal (bounded transition rate) relaxation dynamics 
cannot keep pace with the rate energy is added to the 
system. While an actual divergence of x may be im- 
peded by various physical mechanisms (which we did not 
attempt to incorporate), we believe that the existence of 
a change of behavior at some high strain rate is phys- 
ically meaningful, and likely to be an important prop- 
erty that it is closely related to the existence of a la- 
tency time in stop-start tests. Because transformation 
rates are bounded, at high strain rates, the underlying 
state variable is expected to decouple from the relation 
between stress and strain rates. This decoupling and 
has several consequences: (i) since the shear rate is only 
weakly dependent on state variables, and since it is an 
increasing function of the stress, velocity strenghening 
behavior emerges, (ii) The state variable is hidden, its 
value cannot be deduced from the relation between stress 
and strain rate, and measurements of apparent steady 
relationships between stress and strain rate may not cor- 
respond to a true steady state, (iii) The hidden value 
of the state variable depends on the strain, (iv) When 
the deformation is interrupted, the state variable requires 
some time before it reaches a sufficiently low value for it 
to make a measurable impact that can be probed by the 
emergence of a transient peak (i.e. a stiction spike) in 
stress. In this case, the latency time becomes a func- 
tion of the overall strain. If a high strain rate fixed 
point is reached by the state variable, this dependence 
disappears, and the latency time should only depend on 
strain rate. It would therefore be particularly interesting 
to see whether the underlying dynamics of a state vari- 
able can be probed in experiments close to and above e* 
and whether it is consistent with either the emergence 
of a high strain rate fixed point, or with unstationary 
dynamics of the state of the lubricated contact. Veloc- 
ity strengthening has not been observed with the SFA, 
probably because of intrinsic limits due to the finite scan 
length. Therefore, it is likely that the high strain rate 
regimes we discuss here cannot be directly accessed with 
this experimental setup. However, velocity strengthen- 
ing has been observed in numerics, Q and recently in 



the experiments |3J| by Bureau and coworkers. This lat- 
ter experimental set-up is constructed to allow sliding 
over longer distances, and provides access to higher drive 
velocities than those available with the SFA. It thus may 
present opportunities to probe some phenomena which 
arise at high velocities. 

Finally, we note that the constitutive equations pre- 
sented here are derived from heuristic assumptions, which 
extend STZ theory to the rearrangements occurring in- 
side the sheared material in a lubricated contact. Al- 
though our approach is phenomenological, it relies on a 
specific picture of the microscopic and mesoscopic mech- 
anisms of deformation, and thus forms a bridge between 
macroscopic, empirical, descriptions and a fundamen- 
tal understanding of the microscopic physics. The STZ 
picture of plastic rearrangements occurring in localized 
zones is supported by numerical observations [fifij . How- 
ever, it would be useful to also image rearrangements 
in experimental situations. Although the SFA provides a 
very well-controlled environment at the microscopic level, 
it seems unlikely that such atomistic imaging of the lu- 
bricant could be performed. Other materials are more 
likely candidates for such observations. Indeed, the pic- 
ture of elementary rearrangements which emerges here 
is not limited to lubricants, but is expected to apply to 
a wide range of amorphous materials, in particular col- 
loidal suspensions and granular materials. Indeed, recent 
experiments by Gollub have imaged local rearrangements 
in granular bead packs forming the interfacial material 
in friction measurements. These experiments also ex- 
hibit stick-slip at low drive velocities, and a transition to 
steady sliding as the drave rate is increased. A 
further theoretical challenge associated with these mea- 
surements arises because slip does not occur homoge- 
neously in the material, but rather is primarily restricted 
to a relatively narrow dilatent region near the top of the 
bead pack. We expect the results presented here may 
extend to other types of amorphous materials, but in sit- 
uations when strain localization does not preempt the 
application of a homogeneous description. Colloids and 
granular materials may provide simultaneous access to 
macrocopic rheological properties, and microscopic imag- 
ing of elementary rearrangements. 
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dation, and the NSF Grant No. DMR-9813752, and 
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plex Networks. 



APPENDIX I: HOPF ANALYSIS 

We present here the details of the calculation of the 
Hopf bifurcation. For the stress-linear free-volume con- 
stitutive equations (|14f> and Q15fl, The complete dynam- 
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ical studied can be rewritten as 

1' 



X 



A* (e — Eq> exp 
— Ei exp 



X 



a Eq exp 



(35) 
(36) 



The calculation of the Hopf bifurcation point is straight- 
forward: it amounts to considering the trace of the Jaco- 
bian of this dynamical system. However, the calculation 
identifying the critical Hopf point p, C rit dividing the sub- 
and super-critical bifurcation lines requires more lengthy 
calculations and benefits from the introduction of simpli- 
fying notation. Let us write the equations as: 



X = -fzix) +"/i(x)c 2 



(37) 
(38) 



and perform the analysis in this more general framework. 
The stationary solution is determined by 



and 



v 2 = S{ X ) = 



E( X ) 



h(x) 

otfxix) 

h(x)Hx) 



(39) 



(40) 



where the functions E and S have been introduced for fu- 
ture convenience. In the case of equations (|35|) and (|36|l . 
these functions are: 



E(X) 



EqEi 



■ exp 



1 



and 



s(x) = -|~ ex P 

aE 



K — 1 

X 



The Jacobian of this dynamical system reads: 



J = 



-vfi(x) 

2<r/i(x) 



(41) 



(42) 



(43) 



and around the stationary solution, the eigenvalues of the 
Jacobian satisfy: 



A + A 



' E( X ) 

s(x) 



(p + aS>(x)) + a[iE'( X )=0 . (44) 



The Hopf bifurcation occurs when 
M = -aS'(x) 



(45) 



at any point where E(x) is strictly increasing. In our 
case, E is always an increasing function of x, and S is a 
decreasing function of x if and only if k < 1. If S is an 
increasing function of x ( K > 1) _ which also means that 



S o E^ 1 is an increasing function of e, there is no Hopf 
bifurcation, and the steady sliding motion is stable. If 
k < 1, S is a decreasing function of x\ f° r an Y ^ there is 
a critical value of /x where the system undergoes a Hopf 
bifurcation: 

Mhopf = -aS' {E-\e 2 )) 

Ei 1 — k I ae~ 



E (k + 1) 2 \EoEi 



hi 



a e 

Eq Ei 



In order to determine the type of Hopf bifurcation 
(super- or sub-critical), we write our non- linear system 
of ODE's in normal form. For this purpose, we look for 
the linear operator T which transforms the Jacobian as: 



T- 1 JT = 



P 



(46) 



where the eigenvalues of the Jacobian are A± = p ± ioj. 
We also have 



! E(x) 

s(x) 



(47) 



and 



co 2 = -p 2 + a f iE , ( X ) ■ (48) 
The complex eigenvectors associated with A± are 



u± 



1 



(49) 



and T is obtained from the real and imaginary parts of 
these complex eigenvectors, u r = Re(u+),Ui = Im(u + ) 

T = {u r Ui) 



V 



1 

2 a 



'S(X) 



2a 





and 



T" 



-2a 



E(x) 
Six) 



(50) 



The condition p — determines the Hopf bifurcation, 
and at this point, the transformations T and T _1 reduce 
to, 



T = 



s'(x) 

1 



S'(x)E'( X ) 
E(x) 



(51) 



and 



T = 



EM 



S'(x)E'( X ) 



S'(X)E( X ) 
S(x)E'(x) 



(52) 
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Next, we implement the linear change of variables, 



-l & — &o 



X- Xo 

which leads to the system of ODE's, 



-ui 



+ 



f{x,y) 
g{x,y) 



(53) 



(54) 



The stability coefficient a determines whether the bifur- 
cation is super- or sub-critical. In normal form, this co- 
efficient can be directly obtained from the derivatives of 
the functions / and g: 

a Yg if xxx ~\~ fxyy T 9xxy ~\~ 9yyy) 
+ -j^g ~ (fxy i f xx + fyy) 

~9xy {9xx + 9yy) ~ fxx9xx + fyy9yy^j 
We then obtain, at the Hopf bifurcation point, 
a^E(x)(S'(x) 2 -2S(x)S"(x)) 



J 'xx 
fxy 

fvv 



2S( X )3/2 



mx) 

aS'(x)E>( X ) 



2^fW)VEtx) 
ay^SW) {E'{xY -2E(x)E"( X )) 



9xy 



fxxx = 



fxyy 
9yyy 
9xxy 

whence, 



2E{x)^/Wix) 

a^E(x)S'(x) 2 
S( X ) 3/2 
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We thus obtain a quite simple condition for the vanishing 
of a, which determines the points where the super- or 
sub-critical character of the Hopf bifurcation change: 

S^(x) 1 (MM mx) , 2E"( X ) \ 
S"( X ) 2 { S( X ) E( X ) + E'(x) ) 1 ) 

With the specific functions E and S of our current inter- 
est, the parameter a reads, 



1 



16 X < 



■ exp 



£?i(1-k) (k - 1 - 2k X + 2 X 2 ) 



which admits a single positive zero for k < 1: 
Xcrit = 7j (« + v / 2-2k + k 2 ^) 



(56) 



(57) 



This value defines a critical value for the parameter /i at 
the Hopf bifurcation, 



Merit 
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and the critical strain rate 
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